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Abstract. Aims. The ESO SN la Progenitor Survey (SPY) aims at finding merging double degenerate binaries as candidates 
for supernova type la (SN la) explosions. A white dwarf merger has also been suggested to explain the formation of rare types 
of stars like R CrB, extreme helium or He sdO stars. Here we present the hot subdwarf B binary GD 687, which will merge in 
, less than a Hubble time. 

Methods. The orbital parameters of the close binary have been determined from time resolved spectroscopy. Since GD 687 is a 
single-lined binary, the spectra contain only information about the subdwarf primary and its orbit. From high resolution spectra 
the projected rotational velocity was derived. Assuming orbital synchronisation, the inclination of the system and the mass of 
the unseen companion were constrained. 



Results. The derived inclination is i = 39.3+H °. The mass M 2 = 0.71+j?- 2 ? M Q indicates that the companion must be a white 



-5.6 • ™ m *- -0.21 

dwarf, most likely of C/O composition. This is only the fourth case that an sdB companion has been proven to be a white dwarf 
unambiguously. Its mass is somewhat larger than the average white dwarf mass, but may be as high as 0.93 M Q in which case 
the total mass of the system comes close to the Chandrasekhar limit. 

Conclusions. GD 687 will evolve into a double degenerate system and merge to form a rare supermassive white dwarf with a 
^ ' mass in excess of solar. A death in a sub-Chandrasekhar supernova is also conceivable. 

m ; 
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£^ , 1. Introduction SN la play a key role in the study of cosmic evolution. 
■ They are utilised as standard candles for determining the cos- 
Double degenerate (DD) binaries consisting of two white dwarf mological parameters (e.g. Riess et al. [1998l Leibundgut|200TJ 
(WD) stars experience a shrinkage of their orbits caused by Perlmutter et al. [19991 . There is general consensus that the 
gravitational wave radiation. Sufficiently close binaries with or- thermonuclear explosion of a white dwarf of Chandrasekhar 
^ . bital periods of less than half a day will eventually merge within mass cause s a SN la. The DD merger (Iben & Tutukov [19841 
H I less than a Hubble time. The outcome of such a merger can be a [ s one f two main scenarios to feed a white dwarf to the 
. P. , single compact object like a white dwarf or a neutron star. But Chandrasekhar mass. Progenitor candidates for the DD see- 
the merger event may also lead to the ignition of nuclear burn- nar i have to merge in less than a Hubble time and the total 
ing and trigger the formati on of r are objects such as extre me he- binary mass must exceed the Chandrasekhar limit, 
lium stars (Saio & Jeffery 2002), R CrB stars (WebbinktHHB In genend> DD candidates ^ assumed t0 consist of 
or He sdO stars (Heber 2008; Napiwotzki 2008). The merger tWQ whke dwarfs The discovery of the sdB+W D binary 
of two sufficiently massive C/O white dwarfs may lead to a ^ 1930+2752, however, has highlighted the importance of 
Supernova of type la (SN la, Webbink,[1984l>. such systems as SN Ia progenitor candidates (Bi lleres et al. 

[20001 Maxted et al. l200Tjl Geier et al. [20071 . Hot subdwarf 

~~ ~~~Z ~ „ . stars (sdBs) are core helium-burning stars situated at the hot 

Send offprint requests to: S. Geier, , , . , , ... , , TT , 

., . „ . , , end of the horizontal branch with masses around 0.5 M (Heber 

e-mail: geier@sternwarte.uni-erlangen.de . . wv 

* Based on observations at the Paranal Observatory of the European I™- Although the formation of these objects is still under de- 
Southern Observatory for programme No. 165.H-0588(A). Based on bate > Jt ls general consensus that they must loose most of their 
observations at the La Silla Observatory of the European Southern envelope at the tip of the RGB and evolve directly from the 
Observatory for programmes No. 072.D-05 10(B), 079.D-0288(A), EHB to the WD cooling tracks without ascending the AGB (see 
080.D-0685(A) and 084.D-0348(A). Heber[2009 for a review). 
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Since a high fraction of sdBs resides in close binaries with 
unseen companions (Maxted et. al 2001 ; Napiwotzki et al. 
2004), such systems can qualify as SN la progenitors, if the 
companion is a white dwarf and the lifetime of the sdB is 
shorter than the merging time of the binary. This is the case for 
the sdB+WD binary KPD 1930+2752, which is the best known 
DD progenitor candidate for SN la (Maxted et al. 120001 Geier 
et al. 120071 

Systematic radial velocity (RV) searches for DDs have been 
undertaken (e.g. Napiwotzki 2003 and references therein). The 
largest of these projects was the ESO SN la Progenitor Survey 
(SPY). More than 1000 WDs were checked for RV-variations 
(Napiwotzki et al. l200ni2003l . SPY detected ~ 100 new DDs 
(only 18 were known before). One of them may fulfil the crite- 
ria for SN la progenitor candidates (Napiwotzki et al. 2007). 

This sample includes about two dozen radial velocity vari- 
able sdB stars with invisible companions. As the sdB stars are 
intrinsically brighter than the white dwarfs, the nature of the 
companion is not so easily revealed. 

Most sdBs in close binary systems are single-lined and 
no features of the companions are visible in their spectra. In 
this case the nature of the unseen companion is hard to con- 
strain since only lower limits can be derived for the compan- 
ion masses. Given the fact that sdBs are subluminous stars, 
main sequence companions earlier than K-type can easily be 
excluded. But if the derived lower limit for the companion mass 
is lower than about half a solar mass it is in general not possi- 
ble to distinguish between a main sequence companion and a 
compact object like a white dwarf. Despite the fact that the cat- 
alogue of Ritter & Kolb (2009) lists more than 80 sdB binaries, 
the nature of their companions can only be constrained in spe- 
cial cases (e.g. For et al. 120 101 Geier et al. 120081 1. in particular 
from photometric variations due to eclipses, reflection effects 
or ellipsoidal deformations. 

Results for eight DD systems discovered in the SPY survey 
have been presented in papers I-IV (Napiwotzki et al. 120011 
Napiwotzki et al. 120021 Karl et al. [2003~l Nelemans et al. 120051 
Here we study the sdB binary GD 687, a DD progenitor system 
which will merge in less than a Hubble time. 

2. Binary parameters 

2.1. Observations and radial velocity curve 

GD687 (a20oo = 01 h 10 m 18?5, <5 2 ooo = -34°00'26") was in- 
cluded in the SPY survey, because it was misclassified as DA3 
white dwarf (McCook & Sion ll999l) . GD687 was observed 
twice in the course of the SPY project with the high resolu- 
tion echelle spectrograph UVES at the ESO VLT. The spectra 
have a resolution of R « 37 000 and cover 3200 - 6650 A with 
two small gaps at 4580 A and 5640 A. 

Follow-up medium resolution spectra were taken with the 
EMMI (R * 3400,^ = 3880 - 4380 A) and the EFOSC2 
(R » 2200, A = 4450 - 5110 A) spectrographs mounted at the 
ESO NTT. Reduction was done with the ESO-MIDAS pack- 
age. The radial velocities (RV) were measured by fitting a set 
of mathematical functions (Gaussians, Lorentzians and poly- 
nomials) to the hydrogen Balmer lines using the FITSB2 rou- 




0.4 0.6 
orbital period [d] 



Fig. 1. In this power spectrum - \ogx 2 of the best sine fit is 
plotted against the orbital periods. The region around the best 
solution is shown in the inlet. Confidence limits (lcr,3<x,6cr) 
are marked with horizontal lines (dotted, dashed, solid). 



tine (Napiwotzki et al. 2004). The measured RVs are given in 
TableQ] 

Sine curves were fitted to the RV data points in fine steps 
over a range of test periods. For each period the x 1 of the best 
fitting sine curve was determined. The result is similar to a 
power spectrum with the lowest x 1 indicating the most likely 
period (see Fig. [TJ. The plot shows two very closely spaced 
peaks. We performed a Monte Carlo simulation for the most 
likely periods. For each iteration a randomised set of RVs was 
drawn from Gaussian distributions with central value and width 
corresponding to the RV measurements and the analysis re- 
peated. 

The probability that the solution with the lowest^ 2 and P — 
0.37765 d is the correct one is estimated to be 74%. The second 
best alias period (P = 0.37828 d) has Ay 1 = 9 and a probability 
of 26% to be the correct one (see Fig. Q] inlet). None of the 
10 000 Monte Carlo iterations indicated a period different from 
the two discussed above. The probability that any other solution 
is the correct one is therefore less than 0.01%. Our analysis is 
based on the most likely solution. 

In order to derive most conservative errors for the RV semi- 
amplitude K and the system velocity y we fixed the most likely 
period, created new RV datasets with a bootstrapping algorithm 
and calculated the orbital solution in each case. The standard 
deviation of these results was adopted as error estimate and is 
about three times higher than the lcr-error. The phase folded 
RV curve is shown in Fig. [2] Parameters are given in Table [ffl 



1 Adopting the second best solution only leads to small changes in 
the derived parameters (A* = 124.2±3.0kms~ 1 , y = 37.9±3.2kms _1 ) 
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Table 1. Radial velocities of GD 687 




phase 

Fig. 2. Radial velocities of the primary subdwarf plotted against 
orbital phase. The residuals are plotted below. Filled triangles 
mark RVs measured from UVES spectra, filled diamonds mark 
RVs obtained from EMMI spectra and filled circles RVs mea- 
sured from EFOSC2 spectra. 



2.2. Gravity and projected rotational velocity 

The atmospheric parameters of GD 687 were measured by fit- 
ting metal-line blanketed LTE models to the UVES spectra and 
are given in Lisker et al. (2005). 

In order to derive v rot sin i and the elemental abundances, 
we coadded the observed high resolution spectra after shift- 
ing them to rest wavelength. The spectrum was then compared 
with rotationally broadened, synthetic line profiles calculated 
with the LINFOR program (developed by Holweger, Steffen 
and Steenbock at Kiel university, modified by Lemke[T997]l. 

Due to the wide slit used (2.1") for the SPY survey (see 
Napiwotzki et al. 1200 II for details) the resolution of the UVES 
spectra is seeing dependent in most cases. In order to measure 
v rot sin i accurately, the instrumental profile has to be taken into 
account. We used the ESO archive to obtain the seeing condi- 
tions during the exposure times of the UVES spectra (DIMM 
seeing monitor, Sarazin & Roddier 1990) 



took the average 
(1.13 arcsec) and folded our models with the instrumental pro- 
file. 

The projected rotational velocity was measured simulta- 
neously with the elemental abundances to v 10t sin/ = 21.2 + 
2.0km s _1 using the six strongest metal lines. In addition, an 
abundance and rotational broadening fit was done using two 
prominent helium lines (Hei 4472, 5876) and LTE model 
spectra (Heber et al. 2000). The resulting v ro tsin/ = 19.8 + 
2.2kms~ 1 turned out to be perfectly consistent with the value 



mid-HJD 


RV [kms- 1 ] 


Instrument 


2451737.83523 
2451740.76555 


162.9 ± 4.0 
63.1+4.0 


UVES 


2453337.67383 


-40.5 


+ 5.3 


EMMI 


2453338.59413 


17.4 


±4.1 




2453339.61394 


-87.6 


±4.0 




2453339.67229 


-64.7 


±4.1 




2453340.54050 


130.1 


±4.7 




2454252.87242 


-13.5 ± 


11.0 




2454252.92242 


82.4 


±9.8 




2454254.86756 


136.3 ± 


11.1 




2454476.54549 


162.1 


±6.3 




2454477.56241 


-18.6 


±5.8 




2454477.65241 


136.1 


± 8.6 




2455144.65365 


129.9 


±7.0 


EFOSC2 


2455144.67058 


95.6 


±4.0 




2455144.70520 


25.9 


±9.1 




2455144.72910 


-37.1 


±7.0 




2455144.73909 


-31.9 


±5.1 




2455144.77512 


-74.0 


±7.4 




2455145.53816 


-77.6 


±7.1 




2455145.54772 


-84.4 


±5.4 




2455145.55739 


-83.7 


±4.0 




2455145.56266 


-76.0 


±4.5 




2455145.56792 


-73.4 


±3.6 




2455145.57774 


-63.9 


±5.9 




2455145.58299 


-65.8 


± 8.0 




2455145.58826 


-51.1 


±6.4 




2455145.59991 


-52.4 


±5.3 





and has no significant impact on the qualitative discussion of our re- 
sults. 



measured from the metal lines. Since metal lines are more 
sensitive to rotational broadening we adopt the value derived 
from these lines for our analysis. The parameters are given in 
Tabled 

Our results depend very much on the accuracy of the 
v ro t sin i and logg measurements. A thorough discussion of the 
error in surface gravity is given by Lisker et al. (2005). To 
quantify the v rot sin i error, we carried out numerical simula- 
tions. For this, synthetic spectra with fixed rotational broaden- 
ing were computed and convolved with the instrumental pro- 
file. Random noise was added to mimic the observed spectra. 
The rotational broadening was measured in the way described 
above using a grid of synthetic spectra for various rotational 
and instrumental broadenings as well as S/N levels. Variations 
in the instrumental profile and the noise level were the domi- 
nant error sources. Behr (2003) used a similar method to mea- 
sure the low v rot sin i of Blue Horizontal Branch stars from high 
resolution spectra. The errors given in that work are of the same 
order as the one given here for GD 687 (2.0 km s -1 ). 



2.3. Analysis 

The analysis strategy used is described only briefly. For de- 
tails we refer the reader to Geier et al. (2007, 2008 ). Since the 
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Fig. 3. r e ff-log ^-diagram. The location of the EHB and the evo- 
lutionary tracks (dashed curves) are taken from Dorman et al. 
<fl993l . 

spectrum of GD 687 is single-lined, it contains no information 
about the orbital motion of the companion, and thus only the 
mass function f m = M 3 omp sin 3 i/(M comp + M sdB ) 2 = PK 3 /2nG 
can be calculated. Although the RV semi-amplitude K and the 
period P are determined by the RV curve, M s( jb, M comp and 
sin 3 i remain free parameters. 

Nevertheless, the masses can be constrained by assum- 
ing tidal synchronisation (see also Napiwotzki et al. I2001I I. 
Combining the orbital parameters with an estimate of the sdB 
mass and with the determination of its v rot sin i and surface 
gravity, allows the mass of the invisible companion to be con- 
strained. The mass of the sdB primary is taken from the pop- 
ulation synthesis models (Han et al. 2002 2003 ) that predict 
a mass range of M sc jb =0.37 - 0.48 M Q for sdBs in binaries, 
which experienced a common envelope ejection. The mass dis- 
tribution shows a sharp peak at a mass of about 0.47 M (see 
Fig. 22 of Han et al. 125031 ranging from 0.43 to 0.47 M Q . 
This theoretical mass distribution is consistent with analyses 
of close binary systems (e.g. Geier et al. 2007) as well as aster- 
oseismic analyses of pulsating sdBs (see Charpinet et al. 2008 
and references therein). 

If the rotational period of the sdB primary is synchronised 
the rotational velocity v rot = 2nR s ^/P can be calculated. The 
radius of the primary is given by the mass radius relation 
R - ^M s ^G/g. The measurement of the projected rotational 
velocity v rot sin i therefore allows us to constrain the inclina- 
tion angle i. For the most likely sdB mass M^b = 0.47 M the 
mass function can be solved, and both the inclination angle and 
the companion mass can be derived. The errors are calculated 
by chosing the most extreme values for the input parameters 
within their respective error limits. In order to account for the 
theoretical uncertainity in sdB mass, we adopted the predicted 



mass range for the sdB (0.43 - 0.47 M ) and calculated the 
lower limit for the companion mass under the assumption that 
M;dB = 0.43 M . The error budget is dominated by the uncer- 
tainties in the v lot sin i and logg measurements. 

Table 2. Orbital and atmospheric parameters of GD687. 
fTaken from Lisker et al. (2005). $The merging time is cal- 
culated using the formula given in Ergma et al. d20011 >. 



Orbital parameters 


T [HJD] 


2455144.515 +0.002 


P 


0.37765 + 0.00002 d 


7 


32.3 + 3.0 km s -1 


K 


118.3 ± 3.4kms-' 


f(M) 


0.065 ± 0.006 M Q 


Atmospheric parameters 


7" e ir 


24350 ± 360 Kt 


logg 


5.32 ± 0.07 Kf 


V ro , sin i 


21.2 ±2.0 km s" 1 


Derived binary parameters 



Mi (adopted) 


0.47 M s 


i 


0.25 + 0.02 R Q 

on t+6.2 o 


M 2 


0.71 1| M G 


Mi + Mi 


1.18*}* M e 


►merger 


11.1 x lO'yr* 



3. Nature of the unseen companion 

The mass function provides a lower limit to the mass of the in- 
visible companion of 0.35 M . In the case of a white dwarf 
primary it is impossible to hide the contribution of a main 
sequence star even of the lowest mass in optical/NIR spec- 
tra since these are intrinsically faint. This is not the case for 
sdB stars. A main sequence companion with a mass lower than 
0.45 M can not be excluded because its luminosity would be 
too low to be detectable in the spectra (Lisker et al. 2005). 
This is the reason why the companions' nature still remains 
unknown for most of the x;80 sdB systems in the catalogue of 
Ritter & Kolb (2009). Additional information is needed. 

No spectral features of a cool main sequence star are 
present in the optical spectra of GD687. Furthermore, Farihi 
et al. (120051 1 included GD 687 in a near-infrared imaging sur- 
vey to search for low-luminosity companions to white dwarfs 
and found no evidence for an infrared excess which could be 
caused by a cool main sequence companion. 

As the lower limit for the mass of GD 687's companion de- 
rived from the mass function is lower than 0.45 M , additional 
information is needed to clarify its nature. We made use of the 
gravity and projected rotational velocity to constrain the mass 
to 0.7 1 +y 2i M . The companion therefore can not be a main se- 
quence star but has to be a white dwarf. Taking into account the 
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possible mass range it is very likely to be of C/O composition. 
Its mass exceeds that of an average white dwarf. 

GD687 is only the fourth sdB star, for which the white 
dwarf nature of the companion could be shown unamibigu- 
ously. The others have been discovered by analysing ellipsoidal 
light variations (KPD 1930+2752, PG0101+039) or eclipses 
(KPD 0422+5421, Orosz & Wade fl999l . KPD 1930+2752 and 
PG 0101 +039 have been confirmed with the method used here 
(Geier et al. 120071 120081 . 

4. Discussion 

The derived companion mass was calculated under the assump- 
tion of orbital synchronisation. Since theoretical synchronisa- 
tion timescales for hot stars with radiative envelopes are not 
consistent (Zahn [1977l Tassoul & Tassoul 1992 ), empirical ev- 
idence for orbital sychronisation in sdB binaries is needed. 
Geier et al. (|2008) found such evidence by detecting a vari- 
ation in the lightcurve of the sdB+WD binary PG 0101 +039, 
which could be identified as ellipsoidal deformation of the sdB. 
Since the orbital period of PG 0101 +039 is 0.57 d, sdB bina- 
ries with shorter periods like GD 687 are very likely synchro- 
nised as well. Recently van Grootel et al. (120081 ) performed 
an asteroseismic analysis of the pulsating sdB binary Feige 48 
and for the first time proved orbital sychronisation in this way. 
The orbital period of Feige 48 (0.36 d) is very similar to the 
one of GD687. Furthermore, the atmospheric parameters of 
GD687 (T eff = 24 300 K, logg = 5.32) indicate that it has 
already evolved away from the ZAEHB (Zero Age Extreme 
Horizontal Branch) and should therefore, according to evolu- 
tionary calculations (e.g. Dorman et al. 119931 ), have a similar 
age as PG 0101 +039 and Feige 48 (see Fig. [3). We thus con- 
clude that the assumption of orbital sychronisation is fully jus- 
tified in the case of GD 687. 

In about lOOMyr the helium-burning in the core of the 
sdB will come to an end. After a short period of helium-shell- 
burning this star will eventually become a white dwarf consist- 
ing of C and O. GD 687 is one of only a few known DD pro- 
genitor systems, where both components are C/O white dwarfs 
and which will merge in less than a Hubble time. 

Compared to the sdB + WD binary KPD 1930+2752 (P * 
0.1 d, Maxted et al. 120001 Geier et al. 12007] ), which is the 
best known candidate for DD SN la progenitor, the orbital pe- 
riod of GD 687 is rather long. This leads to a merging time of 
11.1 Gyr, which is just a little shorter than the Hubble time, 
compared to only 200 Myr for KPD 1930+2752. With a total 
mass of 1.18^Q 2j M Q for the most likely subdwarf mass it may 
come close to the Chandrasekhar limit of 1 .4 M Q and is there- 
fore placed at the edge of the progenitor parameter space (see 
Fig-El). In contrast to KPD 1930+2752, where the primary mass 
could be constrained by an additional analysis of the subdwarfs 
ellipsoidal deformation visible as variation in its lightcurve, no 
such constraint can be put on the primary mass of GD 687 yet. 

Instead of exploding as SN la, the merger of the two white 
dwarfs will most likely lead to the formation of a supermassive 
white dwarf with O/Ne/Mg-core. 

Up to now four binaries with total masses between about 
1 .2 and 1 .4 M have been discovered. In two of them the visible 
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Fig. 4. Total mass plotted against logarithmic period of double 
degenerate systems from the SPY survey. GD687 is marked 
with the filled circle. For KPD 1930+2752 a mass range is 
given (Geier et al. 2007). The filled rectangles mark double- 
lined WDs, for which the absolute masses can be derived. The 
filled circle with arrow marks the lower mass limit derived for 
HE 1047-0436 (Napiwotzki et al. BDOTT) . the filled diamond the 
lower mass limit derived for PNG 135.9+55.9 (Napiwotzki et 
al. ©05 ). The open symbols mark single-lined WDs, sdBs, and 
sdOs. The companion masses of the single-lined systems are 
derived for the expected average inclination angle (i = 52°) 
(Napiwotzki et al. 120021 Karl et al. 120031 Karl l2004l Nelemans 
et al. 120051 Napiwotzki et al. 120071 . 



component is an sdB. Since the sdB close binary fraction is 
much higher than the one of white dwarfs, it may be easier to 
find double-degenerate binary progenitors in the hot subdwarf 
population. 

Two other massive DD systems, a central star of a planetary 
nebula (Tovmassian et al. 2004; Napiwotzki et al. 2005) and a 
white dwarf from the SDSS survey (Badenes et al. 2009, Marsh 
et al. 120101 Kulkarni & van Kerkwijk 2010), were discovered 
serendipitously. 

Even though GD 687 does not qualify as SN la progenitor 
candidate, the discovery of a population of double degenerate 
binaries (and progenitor systems) with total masses close to the 
Chandrasekhar limit in the course of the SPY survey (see Fig. 
|4]) provides evidence for a similar population exceeding this 
limit. The same binary evolution channel that produces the sub- 
Chandrasekhar systems will also produce super-Chandrasekhar 
systems with slight changes in the initial conditions only. 

Acknowledgements. S. G. is supported by the Deutsche Forschungs- 
gemeinschaft under grant Hel354/40-3. 
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